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Abstract 

We investigate the existence of invariant measures for self-stabilizing 
diffusions. These stochastic processes represent roughly the behavior of 
some Brownian particle moving in a double-well landscape and attracted 
by its own law. This specific self-interaction leads to nonlinear stochastic 
differential equations and permits to point out singular phenomenons like 
non uniqueness of associated stationary measures. The existence of several 
invariant measures is essentially based on the non convex environment and 
requires generalized Laplace's method approximations. 
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1 Introduction 

The aim of this paper is to present some new and surprising results concern- 
ing the existence of invariant probability measures for one-dimensional self- 
stabiHzing diffusions. The specificity of such diffusion is the attraction of its 
paths by the own law of the stochastic process. The dynamical system solved by 
self-stabilizing diffusions can be characterized by three essential elements: first 
the system is governed by a double- well potential V which represents roughly the 
environment of the process, secondly some interaction potential F describes how 
strong the attraction between the process and its own law is, and finally the sys- 
tem is perturbed by some Brownian motion with small amplitude {^/cBt, t >0). 
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Let us denote by {dx) the law of the self-stabilizing diffusion {X^ , t > 0), then 
the SDE satisfied by (X^) is given by: 

X^ = Xo+y^Bt-f V'{Xl)ds- [ [ F'{Xl-x)dul{x)ds, e > 0. {E'^^°) 
Jq Jo Jr 

Introducing the notation of the convolution product, {E^'^") can be written as 
follows: ^ 

X^=Xo + V~eBt- f {V + F'*ul){XI)ds. (1.1) 
Jo 

Let us just note that the interaction part of the drift term is related to the diffu- 
sion in some simple way: F' *ul{x) = M[F' {x — X^)]. This way of characterizing 
the drift term essentially points out the structure of the attraction between the 
paths of the diffusion and its law. Self-interaction corresponds obviously to 
mean fields stabilization. 

Self-stabilizing diffusion paths can usually be approximated by the movement 
of some specific Brownian particle belonging to a huge ensemble of identical 
ones. In this global system each particle is submitted to the same forces. First 
it moves in the potential landscape characterized by the double-well function 
V and accordingly it is attracted by positions which minimize the potential. 
The second force which acts on the system is the interaction between all the 
particles. More precisely each one is attracted by all the others. This attraction 
can for instance be thought of as being generated by electromagnetic effects. 
In this case, the solution of the global system doesn't represent some spatial 
position but some electromagnetic charge. 

The huge particle system containing TV elements is governed by the following 
stochastic differential equation 

a;o e M, l<i<N. (1.2) 

Here the are independent Brownian motions. In the limit, as N becomes 
large, the interaction part of the drift term is approximatively the average with 
respect to the law of one characteristic particle of the system (law of large 
numbers framework). More precisely the empirical measures '^jfj " con- 

verges to some law for each fixed time and noise intensity, and each individual 
particle's motion converges in probability to the solution of the diffusion equa- 
tion 

dXi = VedW^ - V\Xi) dt- I F'{Xl - x) du;{x) dt. (1.3) 

Jr 

Interacting particle systems such as (|1.2|) have been studied from various points 
of view. A survey about the general setting for interaction (under global Lips- 
chitz and boundedness assumptions) may be found in f^. 



dX, 
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The aim of this paper is to consider both the existence and the uniqueness 
of stationary measures for the self-stabiHzing diffusion {E'^''^°). In [5] Her- 
rmann, Imkeller and Peithmann proved the existence of some unique strong 
solution to equation (jl.ip generaHzing previous results obtained by Benachour, 
Roynette, Talay and Vallois [ij in the context of constant environment poten- 
tial V {V'{x) = for all a; G M). We choose their work as basis for developing 
our study. Nevertheless there exist several different papers dealing with the 
existence problem for self-stabilizing diffusion, each of them concerning other 
families of interaction functions. Let us cite McKean who studied in some 
earlier work a class of Markov processes that contains the solution of the lim- 
iting equation under restrictive global Lipschitz assumptions for the interac- 
tion [6], Stroock and Varadhan who considered some local form of interaction 
[H] , Oelschlager who studied the particular case where interaction is represented 
by the derivative of the Dirac measure at zero [7] and finally Funaki who ad- 
dressed existence and uniqueness for the martingale problem associated with 
self-stabilizing diffusions [3]. 

Let us now focus our attention to the stationary measures. In [1] the authors 
emphasize that the invariant measure, corresponding to some given average, is 
unique in this particular constant potential V situation. This feature is essential 
for further developments. The natural convergence question between the law of 
the process and the invariant measure, as time elapses, can then be analyzed, 
see [2]. This kind of convergence was also considered by Tamura under different 
assumptions on the structure of the interaction, see [11] and |10) . 

The presence of some potential gradient which describes the environment of 
the self-stabiHzing diffusion is essential for the question of existence and unique- 
ness of invariant measures. In particular, if the landscape is represented by some 
symmetric double-well potential then surprising effects appear due to the lack 
of convexity: we shall prove that, under suitable conditions, there exist at least 
three invariant measures of which one is symmetric (Theorem I4.5P and two are 
asymmetric or so-called outlying (Theorem 14. 6p . In the particular linear inter- 
action case {F'{x) = ax with a > 0), these three measures constitute the whole 
set of invariant measures (Theorem I3.2p provided that V" is a convex function. 

The material of this paper is organized as follows: first we list several as- 
sumptions concerning both the interaction function F and the environment po- 
tential V which permit in particular to assure the existence of the self-stabilizing 
diffusion {E'^''^°). In Section 2 preHminary results concerning the structure of 
the invariant measure (if it exists !) are developed. These results are essential 
for the construction of such measures. The question of existence starts to be 
addressed in Section 3 in the particular linear interaction context. After point- 
ing out some symmetric and asymmetric invariant measures, we point out some 
nice context for which the whole set of stationary measures can be described. 
This study is finally extended to the general interaction case in the last section. 
We postpone different tools concerning asymptotic analysis based on Laplace's 
method to the Appendix. 
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1.1 Main assumptions 

In order to study invariant measures for self-stabilizing diffusions, we especially 
need that (|l.ip admits some unique strong solution. For this reason, we assume 
that both the potential landscape V and the interaction function F satisfy 
some growth conditions and some regularity properties. Moreover we add some 
technical assumptions which permit to simplify the statements. 
We assume the following properties for the function V: 

(V-1) Regularity: V G C°°(M,R). C°° denotes 
the Banach space of infinitely bounded con- 
tinuously differentiable function. 

(V-2) Symmetry: V is an even function. 

(V-3) y is a double-well potential. The equation 
V'{x) = admits exactly three solutions : 
a, ~a and with a > 0; V"(a) > and 
V"{0) < 0. The bottoms of the wells are 
reached for x = a and x = —a. 

(V-4) There exist two constants C4, C2 > such 
that Vx G R, V{x) > Cix^ - Cix'^ . 

(V-5) lim V""(a;) = H-oo and Vx > a, V""(a;) > 0. 

X — >±oo 

(V-6) The growth of the potential V is at most polynomial: there exist 9 S N* 
and > such that \V' {x)\ < C, (l x^i). 

(V-7) Initialization: ^(0) 0. 

Typically, F is a double- well polynomial function. But our results can be applied 
to more general functions: regular functions with polynomial growth as |a;| 
becomes large. We introduce the parameter ■& which plays some important role 
in the following: 

I? = sup-y"(a;). (1.4) 

Let us note that the simplest example (most famous in the literature) is V{x) = 
4 2 

\ — ^ which bottoms are localized in —1 and 1 and with parameter d ~\. 
Let us now present the assumptions concerning the attraction function F. 

(F-1) F is an even polynomial function. Indeed we consider some classical 
situation: the attraction between two points x and y only depends on the 
distance F{x — y) = F{y — x). 

(F-2) F is a convex function. 

(F-3) F' is a convex function on ]R+ therefore for any x > and y > such 
that x>ywe get F'{x) - F'{y) > F"(0)(x - y). 
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(F-4) The polynomial growth of the attraction function F is related to the 
growth condition (V-6): \F' (x) - F' {y)\ < C,|x - y|(l + |x|29-2 + |y|29-2). 

Let us define the parameter a > which shall play some essential role in fol- 
lowing: 

F' (x) ^ ax + Fo{x) with a = F"(0) > 0. (1.5) 

In [5], Herrmann, Imkeller and Peithmann present sufficient conditions for the 
SDE l|l.ip to admit a unique strong solution. In particular, if E[Xq'^ ] < +oo, 
with q defined in (V-6), and if all the main assumptions just defined are satisfied, 
the existence and uniqueness of the solutions are proved. In the following we 
will always assume that the {8q^)-th moment of the initial value Xq is finite. 
This permits to study further the self-stabilizing diffusions and exhibit invariant 
measures. 



2 General structure of the invariant measures 



This section deals with different preliminary results describing the main struc- 
ture of the invariant measures of (£"'^«). First of all, there is some classical 
Hnk between the stochastic differential equation and the associated paraboHc 
partial differential equation which permits to characterize stationary measures. 

Lemma 2.1. Let u^{x) denote the density of (X^; t > 0) with respect to the 
Lehesgue measure. Then v!^ is solution of the following PDF: 



(2.1) 



for allt> 0, X £M. and ul{dx) = P {Xq G dx). 
We recall that J^x^'' UQ{dx) < co. 



Proof Let f € 



I) such that 

lim f{x)^ lim f'{x)=0. 

>-±oo X- >±oo 



(2.2) 



By Ito's formula, we obtain 

E[f{Xt)] = E[/(X^)]+E 



fiXDVedB., 



E 



f'{Xl)(v'iXt) + F'*ul{Xt))ds+'- I f"{Xl)ds 



Taking the time derivative, we get 
d 



dt 



E[/(X| 



-E 



nx^)(v'ixt) + F' * utix^)) + '-nxt) 

nx){y'{x)+F' ^ul{x))ul{x)dx~ j ^f"ix)ut{x)da 
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Since / is a function, integration by parts leads to 



HI 



{V'{x) + F' * ul{x)) ul{x)] + \^ul{x) \ dx. 



Using the equality: 



— / f{x)ul{x)dx 



f{^)-Qlu't{x)dx, 



we deduce for all C'^-functions satisfying (|2.2p : 



[ 



\ dx 



{V'{x)+F' *u\{x))ul{x) 



e d'^uKx) 
2 dx^ 



We obtain (|2.ip by identification. 



dx 



□ 



The density of (Xj, t > 0) with respect to the Lebesgue measure is solu- 
tion to the parabolic PDE (|2.ip (non-linear Kolmogorov equation): this impHes 
in particular that any stationary measure (if it exists !) satisfies some eUiptic 
differential equation. This link between non-linear differential equations and 
self-stabilizing diffusions permits to express the invariant measure in some ex- 
ponential form. 

Lemma 2.2. // there exists an invariant measure to {E'^'-^°) whose (Sg^)- 
moment is finite, then: 



Ue{x) 



1 



A(«.) 
1 



exp 
exp 



F' *u,(y)dy + V{x) 
(f *Ue{x)~ F* u,{0) + V{x)^ 



(2.3) 



where A(ue) denotes the normalization factor: J^Ug{x)dx = 1. 
Proof. By (|2.ip . any stationary measure satisfies 

|<(x) -I- (u,{x) {V'ix) + F' *u,{x))'^' = 0, for all xGR. 

By integrating the previous equality, we obtain the existence of some constant 
Ce e M such that 

|<(a;) + u^{x){V'ix) + F' * u,{x)) ^ C„ for all x G R. 

Using the method of variation of parameters, the solution Uc takes the following 
form 

2 



with 



Ue{x) = A(:(x) exp 
2 



K{x) 



■ Ce GXp 



F' * Ue{y)dy + V{x) 
F' *u,{y)dy + V{x) 
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Hence 



Ueix) 



Ac (0) cxp 



F' * u^{y)dy + V{x) 



e 



exp 



exp 



F' * u,{z)dz + V{y) 
F' * u^{y)dy + V{x) 



dy 



Let us assume that Cc ^ 0. Applying Lemma lA.ll to the function U{x) = 
F' * u^{y)dy + V{x), whose second derivative is positive for |a::| large enough 
(using hypotheses (V-5) and (F-2)), permits to exhibit the equivalent of Ae(a;): 



e " 2 1 



Hence 



Ue{x) 



{V'{x)+F' *u,{x) 



as X 



±CX3. 



V'{x)+F' ^u,{x) 

Due to the conditions (V-6) and (F-4), there exists some constant K > Q such 
that 

\V'{x) + F' * u,{x)\ < K{1 + |.t|-*«-1), for all xeR. 



We deduce that : 



can't be integrated with respect to Ue'- that contradicts 



the essential assumption of the statement. We deduce that Ce = and obtain 
II2.3D after normalization. □ 



Lemma 12.21 presents the essential structure of any invariant measure. The 
global exponential form will play a crucial role in next sections: to prove the 
existence of some stationary measure, it is necessary and sufHcient to solve 
equation (|2.3p . 



3 The linear interaction case 

First we shall analyze the existence problem for stationary measures in the 
simple linear case. In this case F'{x) = ax with a > 0, the interaction gradient 
function is quadratic: F{x) = ^ x'^ and the stochastic differential equation takes 
an interesting simple form. The non-linearity of the drift term is limited to the 
average of the density uf (x) : 

= Xa + yf^Bt - [ V'{Xl)ds -a I [xl- [ xdul{x)) ds, e > 0. 
Jo Jo ^ Jm ' 

The study of this particular case emphases the existence of several invariant 
measures. The interesting problem is then to determine in which situations the 
number of such measures is perfectly known. 
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3.1 Existence of invariant measures 



The existence question is really simplified in the linear interaction case, it is 
just reduced in fine to the following parametrization problem. Let us denote 
the first moment of an invariant measure Ue by 



(3.1) 



then (12.31) becomes 



Ueix) 



exp 


-| (v{x) + - ami{e)x^ 






(V{y) + ~ ami{e)y^ 


dy 



(3.2) 



We now come to the essential equivalence: is an invariant measure if and 
only if l|3.ip and l|3.2p are satisfied. It suffices then to point out the convenient 
parameters mi(e) since there is a one to one correspondence between these 
parameters and the invariant measures. In other words, we shall find the solution 
of the equation 



^e{m) with ^e{m) 



J^xexp[~^(vi 



dx 



c)] dx 



(3.3) 



Obviously, mKe) = is a candidate. The corresponding measure u° is invariant 
and symmetric: 



exp 



-I V(x)+4 



exp 



dy 



In fact M° is the unique symmetric stationary measure. 

Of course the natural question concerns the existence of others reals mi(e) solu- 
tions of l|3.2p . In fact the basic dynamical system associated to self-stabilizing 
diffusions is symmetric since F and V are assumed to be even functions. The 
consequence is immediate: if the initial law of the diffusion {X^, t > 0) is 
symmetric so will be the law of for all t > 0. In [l], the authors consider 
self-stabiHzing diffusions without the environment potential V. They proved 
the existence of some unique symmetric invariant measure and describe the be- 
havior of the diffusion: for any initial law satisfying the moment condition of 
order 8q^ the law of Xt — IE[^o] converges to the invariant symmetric law as 
time elapses. 

Adding some double-well potential V in the main structure of the stochastic 
differential equation changes drastically the situation. In particular we prove 
the existence of several invariant measures, one of them being symmetric. 

Proposition 3.1. Let a be the unique positive real which minimizes V (see (V- 
3)). For all 5 e]0, 1[, there exists eo > such that for all e < cq, the equation 
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I|3.3p admits a solution satisfying the estimates: 



a- . / r^TTf^ e< mi(e) < a- — ^— )-f-- e. (3.4) 



4F"(a) (a + F"(a)) 
Moreover -~mi{e) satisfies (|3.3p too. 



W'ia) {a + V"{a)) 



Let us note that, for e small enough, the preceding proposition impHes the 
existence of at least three invariant measures corresponding to the averages: 0, 
mi(e) and — mi(e). 

Proof. Set r > 0. Let's proceed to the first order asymptotic development of 
the expression ^'c(a — re). 



a{a — Te)x 



dx 



Jjj cxp [— I (T^(a;) + — a(a — re)x)] dx 



xe ^""^"^ exp 




V{x) + — aax 




dx 


4e-2"-exp 


— 1 (F(x) + — aax)] dx 



By Lemma Fa. 51 applied to the context: f{x) = —2aTX, n = 1, U(x) = V(x) 
^x^ — ax and /i = 0, we get: 

1 



4'e(a - re) 



4a {a + V"{a)) 

V"{a) 

re H 

a + V"{a) 



aV^^\a) + Aaar {a + V" (a)) e + o(e) 
F(3)(a) 



4:V"{a) {a + V"{a)) 



e + o(e). 



Then a — r"e is the first order approximation of the 



Oet r - 4V"{a){a+V"{a)) 

fixed point. Indeed for 6 e]0; 1[ we can define 

d± := (a - r°(l ± S)e) -{a- r°(l ± S)e) = ±S ^"j^^^ . r"e + o(e). 

For e small enough, c?+ > and < 0. Since the function is continuous, 
there exists mi (e) G [a — r''(l + (5)e; a — r°(l — (5)e] which satisfies ^'e(mi(e)) = 
mi(e). Finally, by the change of variable x := —x in the integral expression 
IjS^ . we obtain '^e{-mi{e)) = -^eirni{e)) = -mi(e). □ 



3.2 Description of the set of invariant measures 

According to Proposition 13.11 we know there are at least three invariants mea- 
sures. One of them is symmetric corresponding to the average and two others 
will be called outlying measures, one wrapped axound a and the other one around 
—a. The aim of this section is to study if there are exactly three invariants mea- 
sures or more. 

For this purpose, we study the asymptotic behavior of the function ^'e defined 
by l|3.3|) in the small noise limit. 
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Theorem 3.2. // V" is a convex function then, in the small noise limit, there 
exist exactly three stationary measures. 

Proof. Let m > 0. Let us recall that the interaction function is linear: F'{x) = 
ax with a > 0. In order to study the invariant measures, we have to consider 
the fixed points of the apphcation ^'^(to) defined by l|3.3p . We introduce the 
following potential function: 

Wm{x) = V{x) + — — amx. 

Since ^'(0) = 0, we have M^™(0) < 0. Moreover liuix^+oo Wln{x) = +oo. 
So we denote by Xm, the positive real for which the potential Wm, admits its 
global minimum. It is uniquely determined since V" is a convex function. In 
particular, Xm satisfies V'{xm) + ce{xm — ™) = and V"{xm) + a > 0. Fur- 
thermore V"{xm) + a > 0. Indeed, since is a global minimum, the equality 
V"{xm) + a = implies that V'-^^Xm) = that is Xm = which contradicts 
the assumption concerning the positivity of x„i. 
We define 

Xt{m) = ^ ^{m) - m and xo(to) = a^m - m. 



We obtain the expression: 

-^m) exp 

Xe(TO) Xra - m + ■ 



{V[x) + ^ 



dx 



4exp[-|(y(x) + ^ 



c)] dx 



(3.5) 



It suffices to prove that Xe has just one zero in E!j_. 

Step 1: For all e > and m > 0, we observe that Xe("i) < Xo(™) = Xm — rn. 
We apply the change of variable x := y + Xm to the integrals in p.5p and obtain 

Xe{m) = Xo("i) + 
= Xo{m) + 



/rI/cxp 




y{y + Xm) + ^- + a {xm - m) y^ 


dy 


/nexp 


-| (viy + Xm) + ^ + a {xm - m) y^ 


dy 



/o°°2/exp [-fy^] n^^m{y)dy 



/jj exp -| {v{y + Xra) + ^ + a - m) y^ dy 



with 



^e,ni{y) = exp 

— exp 

We introduce the function 



[V{y + ) + a {xra - m) yj 

^ (^(y " ^m) - OL {Xra - m) y 



^m{y) = V{y + Xra) " V {y - Xra) + 2a (x„ - m) y 
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Since V is an even function, Am(0) = and A''„(0) = 0. According to the 
definition ofxm, AJ„(0) = 0. V" is a convex function therefore V^^^^ is increasing. 
So A^^(y) = V(3)(y + Xm) - V^^'>{y - a;™) > for all y. We deduce that A'^ is 
increasing. Hence A"„ is nonnegative on so does Am{y) for y > 0. Finally we 
get Qe,m{y) < for all y > 0. We obtain the announced result: Xe(™) ^ Xoim) 
for m > 0. 

Step 2. xo has a unique zero on W*^. 

Let us compute Xo(a) with a defined in (V-3). We know that a is solution of 
V'{x) +a{x-a) = with V"{x) + a > 0. Hence Xo(a) = 0. 
Let us focus our attention to the variations of the function Xo on the interval 
]0,a]. Since V {x„i) +axm = ctm, and a + V" (xm) > we deduce that m Xm 
is derivable; we obtain 



and 



X" ^™ — TITTTT V > ^ "^^^^^ implies Xo(»7i) ^ }n,T^ v (3-6) 

The denominator is positive due to the definition of Xm- According to (V - 5), 
V"{x) > for all x > a. Hence Xo('^) < for m > a. Since Xo(a) = we 
deduce that, for all m > a, xo("t-) is strictly negative and therefore the function 
Xo has no zero on ]a; +oo[. 

It remains to study xo on the interval ]0, a]. Since V" is a convex function, we 
deduce that the derivative of xo is non positive for Xm > c with c > satisfying 
V"{c) = 0. We know that c > is unique since V"{0) < and V" is a convex 
function. Moreover c < a. Since the function m Xm is increasing for m > 0, 
we deduce that x'o is negative for x g] max(0, mc), a] where rric — c + By 
construction, if rUc > then the equality Xma = c holds. 
We observe then two different cases: 

• If mc < i.e. a < 1^^^--^: xo is decreasing on with xo(a) = 0. The 
unique zero of xo on M*^ is a. 

• If rric > then xo, which is a continuous function on W*_^, is increasing on 
]0,TOc[ and decreasing on ]mc, +oo[ with Xo(a) = 0. It suffices to prove 
that lim„j_>o+ Xo(™) > in order to conclude that a is the unique zero 
of Xo on M*^. Due to the definition of Xm we get: lim,„_+o+ Xo(w) = 
lini,„^o+ > 0. Indeed m Xm is continuous from ]0, +oo[ to ]0, +oo[ 
so the extension to m = is non negative. 

In these two cases, there is a unique zero of xo on M^. 

Step 3. The family of functions (xe)e (respectively (Xe)e^ converges uniformly 

towards xo (resp. 'x!q) on each compact subset o/R^. 

First we prove the convergence of Xe(™) for m > 0. Recall that 



/jjxexp -^(v{x) + ^-amxj 

Xeim) - 



dx 



/gBxp [— I (y{x) + ^1 amx)] dx 
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By Lemma rA.Sl with U{x) — V{x) + n = 1, /i = m and G ~ —ax we obtain 
the announced convergence result: 

Xe{m) - xo{m) = Xeim) - x„, + m = - ^ ^ ^ + 

Ala + V"[Xjn)) 

Moreover this convergence is uniform with respect to the variable m on compact 
subsets of R!j_. 

We estimate now the asymptotics of Xe(™) as e becomes small. Taking the 
derivative of ^'e , we obtain 



5f^(m) = — 
e 



2a I /jgX^exp [-|Wm(x)] dx I J^xexp [-|Wm(a;)] dx\ 



J^cxp[-^W,nix)] dx y /jjCxp da; 

riance of the n 
to the average m by (|3.2p . Hence 



We recognize the variance of the measure Ue™^ which is the measure associated 



X',(m) = ^Var(«('")) - 1. (3.7) 

Applying again Lemma fA.Sl with U — V{x) + G = —ax, ^ = m and n = 2, 
we obtain 

"^^T r^T^ — 7^ ^xi, - ^ ^ — — e + o(e). 

4exp[-fM^„(x)]dx " 2{a + V"{xm)f 

Applying the same lemma with n = 1 permits to compute the first moment: 

/gxexp [-f W^(a:)] dx ^ ^ XmV'^^\x„,) ^ ^ 

By (|3.7p and the computations of the two first moments, we get 

= + = + "(I)- (3.8) 

a ^ V \Xra ) 

Furthermore this convergence is uniform with respect to the variable m on com- 
pact subsets of M*^. 

Step 4. For any S > small enough, there exists eq > such that Xe has a 
unique zero on [S,(X}[ for all e < eq. 

Since there is no zero of Xe on the interval ]a, +oo[ (Step 1 and 2), we focus 
our attention to the interval ]0, a]. On each compact subset of this interval, Xe 
converges uniformly towards the Hmit function xo (Step 3). Hence the zeros 
of Xe are in a small neighborhood of the unique zero of xo namely a (Step 2) . 
Let us study the derivative of Xe in a neighborhood of a. Since x'e converges 
uniformly towards Xo (Step 3) and Xo('^) < in a neighborhood of a (Step 2), 
we obtain that Xc("^) < in a neighborhood of a for e small enough. Finally 
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we proved that, as soon as e is small enough, the function Xe can't admit two 
zeros or more on 

Step 5. There exists S > and eq > such that Xe doesn't vanish on ]0,S] for 
all e < Co- 
in this last step, we have to distinguish three different cases depending on the 
values -d and a defined by l|1.4p and Ijl.Sp . 

Step 5.1. We assume a < 1). In this particular case Wo{x) = V{x) + ax'^ /2 
reaches a unique global minimum on R+ for a; = xq > 0. 

Let us fix some small S > (depending on xq: we shall precise it in the follow- 
ing). We prove that, for e small enough, Xeim) = ^eim) — m > on ]0, (5]. By 
the definition of see p.3p . it suffices to prove that N^{m) > for m sjO, S] 
where 



Ne{m)= I .Tcxp Wm{x) 



dx 



exp 



W^{x) 



dx. 



(3.9) 



Obviously ^^^(0) = 0. Let us prove that is non decreasing. Taking the 
derivative, we get 



mx 



dx. 



This expression is in fact non negative. Indeed, using the symmetry property 
of Wq{x) and the upper bound m < (5, we obtain 



N'^{m) 



2a 



(x- - ±-\ 
\ 2a) 



cosh 



/ 2amx \ 



mx sinh 



/2amx\ 



a 

> - 



Ps{x)e 



'x with Ps{x) ~ X — Sx . 

a 



We split the preceding integral into two parts: the first integral In concerns the 
support [0,2(5] and the second integral I2S the complementary support [2(5, oo[. 
We get N^{m) > f (/o+Isa). 

Since the roots of the polynomial function Ps satisfy 



x± 



d±\IS^ + — ] < 25, 
a ■ 



the polynomial is positive on the interval [26, co [ and can be lower bounded by 
Ps{2S) = 25^ — e/a. Lemma [A. 31 implies the existence of some constant C > 
leading to the following estimate as e 0: 



hs > (26^ - e/a) 



xa + l 



/e e 



-^{Vixo)+axl/2) 



2(5 



(3.10) 



provided that xq > 26 (it suffices then to chose 6 small enough). 

Let us finally focus our attention to the lower bound of the integral term Iq. 

Since the minimum value of Ps is —{S'^/A + e/a) and since W"{0) < 0, we have 



/n > - 



4 



V 4 aJ 



(3.11) 



13 



For S > small enough, V{2S) > V{xo) + axg/2 (since the minimum of V{x) + 
ax'^/2 is only reached for x = Xq). Consequently the negative lower bound of Iq 
l|3.1ip is negligible with respect to the positive lower bound of I2S as e becomes 
small. We deduce that there exists eo such that N'^{m) > for all m € [0,(5] 
and e < eo- Since Nf{0) = we conclude that N^{m) > on ]0, 5] and so is x^. 
Step 5.2. We assume a > d. In this case Wo(x) admits a unique minimum 
reached for a; = and Xm converges continuously to as m ^ 0. Using similar 
arguments as those presented in Step 3, we claim that Xe (resp. Xe) converges 
towards xo (resp. Xo) uniformly on [0,a] as e ^ 0. Due to the regularity of xo 
and by the inequality Xo(0) — — a+v'^o) > we obtain the existence of (5 > 
and eo > such that y^'^irn) > for m € [0, 6] and e < eo- Xe starts in and is 
strictly increasing on [0, S] which implies the announced result. 
Step 5.3. We assume that a = -d. It suffices then to note that x« depends 
continuously on the parameter a. The following results can be directly deduced 
from the preceding case (Step 5.2) by continuity: Xe(0) > and Xe(™) ^ for 
m G [0, S] and e < eo. In fact Xe vanishes for x = and is increasing on [0, S]. 
The inequality Xe(™) > for all m €]0, S] and e < eo is an obvious consequence. 
Conclusion: Step 4 and 5 lead to the existence of eo > such that for all e < eg, 
Xe has exactly three zeros: and two other reals, one in the neighborhood of a, 
the other one near —a. To each of these averages corresponds a unique invariant 
measure obtained by p.2p . □ 

Example: In Theorem 13.21 for all a > 0, as soon as e is small enough, there 
exist exactly three invariant measures. There is a one to one correspondence 
between these measures and their average through p.2p . It suffices to determine 
the averages which are in fact solutions to the equation 

X"(m) := \I'e(m) - to = 0. 

These solutions are really close to the solutions of Xo ("^) = 2:^ — to = in the 
small noise limit. We recall that x^ is the global minimum of 

W^{x) ■.= V{x) + ^x^ ~ amx on K^. 

4 2 

Let us observe these averages in the particular case: V{x) = "I ^ and 

F{x) = fx"^. In this case, we compute the parameter c > which vanishes V" 

and the corresponding parameter TOc = c + We obtain: 

1 3a -2 ^ ^ 

c = —= and TOc = — -1= — . (3.12) 
VS 3v3a 

We shall for this example present graphs of the functions Xo (dotted line) and 
X" for different values of a. We choose e = 1/4. Even if it seems to be not very 
small, this value suffices in this example to observe three invariant measures for 
each interaction parameter value considered. 
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First of all we have to determine the value of which is solution of the 
system (i;"'™): 

X^ + {a- 1)X -am = and SX^ + (a - 1) > 0. 

Its discriminant is equal to 

. , , a^?7i^ (a — 1)^ 

A.M = ^ + ^^. 

We distinguish different cases: 

a = 0: the solution is evident, we get = 1 and Xo(to) = 1 — m for 
TO > and by symmetry Xo(to) = —1 — to for all m < 0. Moreover 

x[S(o) = o. 

a > 1 (Figure [2l): for all to e M, we get Aa{m) > 0. Hence 



The function Xo is C°°-continuous and odd. We observe also that Xq (0) = 
Xo(l) = and Xo'("^c) = with defined by ([3A2)) . Hence Xo is 
increasing on ]0, TOc[ and decreasing on ]toc, oo[. 

• a = 1 (Figure [3]) then Xo ('"^) = ~ The limit function is odd, 
continuous on R and C°° on R*. Moreover the path is increasing for 
TO e]0, TOc[, decreasing for to gJjtic, oo[ with rric = 



• 2 



3 < a < 1 (Figure |4]): the discriminant can be negative. Therefore let us 
define mo{a) such that AQ(TOo(a)) = 0. Then for all m between and 
Too (a), the discriminant is negative and for all to larger than too (a) it is 
positive. We get TOo(a) = 2(1 — a)"^ /{3a^/3). We obtain the following 
function: Xo (0) = and 

ifii^^m) \fme [-TOo(a);0[lJ]0;TOo(Q;)] 
ip2"\m) Vto e] - oo; -TOo(a)] |J[too(q:); +(X)[ 




„ / 1 - a 
2\ cos 




27 



(l-a)3 



TO. 



Let us note that Xo (0^) = V^~~a ^ and Xo (0") = -\/T~~^ ^ 0. The 
function is C°°-continuous on ]0; TOo(a)[U]TOo(a); +oo[ and continuous in 
Too (a). 
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Moreover the function is increasing on the interval ]0, mc[ and decreasing 
for m > rric. The maximum is therefore reached for m = me- We observe 
that TOc < "ig for a g [2/3, 3/4] and rric > TOq for a £ [3/4, 1]. We remark 
also that the increasing part is smaller and the decreasing part is longer 
for smaller values of a. 

Furthermore, the part where Xa is equal to (/sj"'' is longer. 

• a < I (Figure [5]) : the function Xo is defined in the same way as in the 
preceding case. The important difference is that the function is decreasing 
on M*^ since rric defined by (j3.12p is non positive. 




Figure 2: Xo (dotted line) and x? Figure 3: Xo (dotted line) and Xe 
for a > 1 for a = 1 




-1 -0.5 0.5 1 ■ -1 -0.5 0.5 1 



Figure 4: 2/3 < a < 1 Figure 5: a < 2/3 

4 The general interaction case 

We assumed for this study that the self-attraction phenomenon is represented 
by a polynomial function F' , see (F-1). In previous section, we analyzed the 
particular linear situation: F'{x) = ax and proved under suitable conditions 
that there exist exactly three invariant measures in the small noise limit. In this 
section we shall focus our attention to the general case: the polynomial function 
F is of degree n > 2. First we shall present results concerning the symmetric 
invariant measure and secondly we discuss the presence of asymmetric measures. 
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4.1 Symmetric invariant measures 

In the linear case we proved the existence of a unique symmetric invariant 
measure. The result is obvious since it suffices to solve the equation l|3.3p with 
mi(e) = 0. In the general case in order to find the symmetric measure we have 
to solve some equation like p.2p but depending on much more parameters than 
just the mean mi(e). The total number of parameters depends in fact on the 
degree of F. Instead of trying to solve such system, we choose some other kind 
of proof based on a fixed point theorem which permits to prove the existence 
of symmetric invariant measures in even more general cases: the interaction 
function does not need to be polynomial. In [1], Benachour, Roynette, Talay 
and Vallois introduced this method of proof for a self-stabilizing diffusion in 
the constant environment case {V'{x) = 0). This proof can be adapted to our 
situation and is based on the following Schauder's theorem (see for instance [4] 
Corollary 11.2 p. 280): 

Proposition 4.1. Let B a Banach space, C a closed convex subset and A a 
continuous application C — > C such that A(C) is compact. Then A admits a 
fixed point in C. 

In order to use this proposition we introduce some definitions and notations: 

1. Let us choose p > Aq where q is defined in (V-6). 

2. D = {w : R — ■> | v is symmetric and supj.g]j+ (1 + \xY')v{x) < oo}. 

3. B = {/:M — >R; sup^gu (1 + l^^l^) 1/(^)1 < oo}- Let us note that D C B. 
B is equipped with the norm | • |oo where |/|oo = sup^.gj{ (1 + \x\p) |/(a;)|. 

4. For all M > we define the function space Cm as the subset of all non 
negative and even function belonging to B which satisfy: 

/ fix) dx = 1 and sup (1 + \x\p) f{x) < M. 

5. For any function / G D we define the operator: 

AHf)(x) - exp[-f(nx) + /;(f-./)(,)dy)] 



1 

■ cxp 



2 



where Ae(/) is the normaHzation factor. 

6. For any function u e B, we define the moments 7fc(w) = ^T^\x\^u{x)dx 
with < A: < p - 2. 

Let us just point out that Cm is a closed and convex subset of B. Moreover 
we have Ca/ C D C B. The aim of this section will consist in proving that 
the application A'^ is C°(Ca/, CA/)-continuous and that A'(Ca/) is compact. 
Therefore Schauder's theorem implies the existence of a fixed point and as the 
matter of fact the existence of an invariant measure in the function space D. 
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Lemma 4.2. For all u £ Cj\/, we have: 

1- lk{u) < MCl where Ci = 1 + max / r-rd-x for all < k < p — 2. 

0<r<p^2j^ 1 + \x\P 

2. there exists a constant C2 > independent of M such that 

< \ {F' *u){y)dy <C2Mx^{l + x^'') for all x. (4.2) 
2 Jo 

Proof. 1. Let u e Cm then the function x i+\x\p integrable on M since 
k < p — 2. Moreover the definition of <Cm impUes that (1 + \x\p) u{x) < M for 
all X e R. Therefore 

= X YT^d + \^\nuix)dx < M I ^dx < MC.. 

2. Let X > 0. Since u e Cm, u is an even function. By l|1.5p we have 
F'{x) = ax + Fq{x) and (F-3) implies that F' and Fq are non negative odd 
functions so is Fq*u. Using the inequality developed in the statement of Lemma 
4.3 in [1] and the assumption (F-3), we have 

F^(x) < i [F^^ix -y) + F[,ix + y)) for yeM, x>0. 
Therefore, for x > 0: 

nX poo 

{Fl,*u){y)dy = / / {Fl,{y-z) + Fl,{y + z))u[z)dzdy 

^0 ^0 

> / / 2F;,{y)u{z)dzdy > 



Jo Jo 

From the preceding inequality we deduce 

a a 

{F' * u) {y)dy ^ (Fq * u) {y)dy + -x^ > -x^ for all x>0. 
Jo 2 2 

Since Jq {F' * u){y)dy is an even function, we get the inequality for all x e 
3. Due to the symmetry oi F' * u we restrict our study to a; > 0. 



[F' *u){y)dy = \J^ ^ [F'{y-z) + F\y + z)y{z)dzdy. 

According to the assumptions (F-1) and (F-4), F is an even polynomial function 
of degree smaller than 2q with g > 1 . We can therefore write F' as follows 



q-l 

2k+l 



F'{x) = ^ akx 

k=0 
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Therefore defining ^(y, z) = F'{y — z) + F'{y + z) we get 
Hy,z)^yY.a,Y^CllX\y^=z^^--^ 

k=0 j=0 

k 



< y max \ak\2 max } 

0<k<q-l 0<j<q'^ 



y^^z^'^-^^ <Cy(l + y^'^) (l + z^?) 



j=o 

Finally since p > Aq, there exists some constant C > such that: 



J'{y,z)u{z)dz < Cy{l + y^'^) / {l + z'^'^) u{z)dz 



< 



dz 



Cy{i + y'') I ^^((i + z>(z)) 

< C'yM{l + y^i). 
By integration we obtain / {F' *u){y)dy < C2Mx'^{l + x^'') for all x e IR+. 

JO 

Lemma 4.3. There exists Mq > such that for any M > Mq, A^(Ca/) C C 



□ 



Proof. By construction A^u is a non negative even function which satisfies 
J^A'^u{x)dx = 1. It suffices then to prove that: 

sup(l + \x\P)A'u{x) < M. 



By l|4.ip and according to Lemma 14.21 we obtain some lower bound for the 
normalization factor: 



+00 



cxp 



> 



— oo 

+00 



cxp 



V{x)+ {F'*u){y)dy 
{Vix)+C2Mx^ (l + x^^)) 



dx 
dx. 



According to both (V-3) and (V-7), we know that V{x) < for all x G [—a; a]. 
Hence 

--C2Mx^{l + a^'i) 



Kill) > 



cxp 



dx. 



Let us define ^{M) = e^/^{2C2M{l + a^'^))-^/^ then limM-oo C(^) = 0. By 
the change of variable x := £,{M)y and Lemma [A. H the following development 
holds 



exp 



-C2Mx2(l + a29) 



dx = 2f (A/) 



a/e.{M) 



~^'dx 







v o, 



■ exp 
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32C2(l+a^ 



As soon as M is large enough, we have Ae(u) > y/Tr£_{M)j - 

Therefore yJu) — C{e)^/M where C(e) is a positive constant determined by 
parameters of the global system and e. By (|4.ip and the preceding upper bound, 
we prove that 

(1 + \x\P)k''u{x) < C(e)VM(l + |x|P)e-t^(^) < C'{e)VM, 

where C"(e) is a positive constant similar to C(e). In order to conclude, it is 
sufficient to choose M > C"(e)^: we get immediately A'^u S Caj. □ 



Lemma 4.4. 

norm. 



is a continuous operator on Cm with respect to the uniform 



Proof. We shall find some upper bound for the following expression |A'^u — A*^?;! 
Step 1. Let u,v G Cm- We define: 



A^(:,)=e-T^(-) 



exp 



exp 



{F'*u){y)dy 
{F^^u){y)dy 



— exp 
exp 



(F' ^v){if)dv 
{F'^^v){y)dy 



It is well known that |e ° 



< |a — 6| for a,6 > 0. In order to apply 



this inequality we have to prove that Jq{Fq * v)[y)dy and Jq{Fq * u)[y)dy are 
non negative. By Lemma 221 for each function / G Cj\/ the convolution term 
{F' * /) {y)dy is lower bounded by fa;^. So Jg" (Fq * /) (y)dy is non negative 
due to the relation: F'{y) = -Fo(y) + ^V- Hence 



|A^(x)| < -e-?^(^)-T-^A^(x), 



(4.3) 



with Aq defined by 

{F^^u)iy)dy^ 



{F^*v){y)dy 



Kiy " ^)("(^) - v{z))dzdy 



Since w and v are elements of Cm , they are even functions and the integral with 
respect to the variable z becomes 



< 



Jo 

X poo 



{Fq{z + y)- Fq{z - yj) {u{z) - v{z)) dzdy 
\F^iz + y)- F^{z - y)\ \uiz) - viz)\dzdy. 



(4.4) 



The assumption (F-4) gives informations about the increments of the interaction 
function: there exist two positive constants Cq and C such that 



22,-l|y|2,-2-) 



\F(,iz + y) - F[,iz - y)\ < 2\y\a, {l + \z + y]'"-' + \z - y\''^~') 

< 2|y|C, (1 + 229-1 |zp9-2 

< c|j/|(i + M2«-i + |zp9-i) 

< C\y\{l + \y\'^-'){l + \z\'^-' 



(4.5) 
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We shall now find some upper bound for \u{z) — v{z)\ in l|4.4p . Since u,v G Cm 
then u{z){l + \z\p) < M and v{z){l + \z\p) < M,yz G R. The obvious upper 

bound \u{z)-v{z)\{l + \z\P) < 2M permits to obtain y/\u{z) - v{z)\ < Jj^^- 



2M 
1 + 



where 



Consequently, for all z of M, \u{z) — v{z)\ < y^\\u — i;| |c 

denotes the uniform norm. Using this inequality, (|4.5p and (|4.4p in order to 
estimate Ag, we get 



\K{x)\<c^\[^ 



|y|(i + |y| 



dy 



/ 2M 

yi+w 



(l + z^'^-')dz. 



Since p > 45 the integral with respect to the variable z is finite and can be 
considered like a constant term. By (|4.3p and using the positivity of ax^, we 
obtain directly the existence of some positive constant C > such that 

M 



According to (V-2), the expression (l + jxp"^"^) e~?^*^^^ can be bounded by 
some constant independent of e. Therefore 



IIA^IU < C(Af, e)^\\u-v\\^. 

Two results can be deduced: firstly 1 1 A*^ | |oo is finite and secondly 1 1 A*^ 
small as \\u — v\\oo decreases towards 0. 
Step 2. For any x e R, we introduce: 



n^ix) 



1 



■ exp 



(F' ^v){y)dy + V{x] 



Then the difference A'^u{x) — A'^v{x) can be decomposed as follows 

1 



A^u(x) - A'wfx) = , , , 

Taking the uniform norm, we get 

1 



A^(:e) + (A,(w)- Ae(u))a(x). 



becomes 



(4.6) 



(4.7) 



Ae(w) 



||A^||oo + |Ae(z;)-Ae(u)|||a||, 



(4., 



We have shown in the proof of Lemma l43l that — C!{e)VM and moreover 
IIA'lloo < C{M, e)^||w- w||oo. We deduce that 
1 



m\oo <C'iM,e)^\\u-v\U 



It is then sufficient to find a similar inequality for the term |Ae(?;) — Xe{u)\ ||r2e 
in order to conclude the proof. 



\K{V)-Xe{u)\ 



A''{x)dx 



IM 



< C\l — V\\u-v\ 



+CXD 



(l + |a;|2*-i)e-*^(^)dx. 
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According to (V-4), the integral with respect to the variable x is finite and does 
not depend on M. We have immediately 



\\,{v) ^ Kiu)\ < CiM,e)y/\\u- v\\^. 
It remains to estimate fleix). By (V-4) and l|4.2p . we have 

{F' * v) {y)dy + Vix) > dx^ + (f - ^a)^^ 

for all X positive. Furthermore the symmetry property of V and F permits to 
extend the bound to all x e R. The function exp [— | {Jq{F' * v){y)dy + V{x))] 
is then bounded by a constant depending on e. Moreover we have already proved 
that j^j^ < C{e)VM for all elements / of the function space Cm- This bound 
can therefore be applied to u and v. Finally we obtain the existence of some 
constant C(e) > such that, for all real value x, \fle{x)\ < C{e)M. 
By (|4Ji) . we have 



IIA^u- A^IU < C'(M, e) Vlk - f lloo + C{M, e) ^11" - w||ooC(e)M. 

In other words, 

IIA^u- A^lloo < C"iM, e)v/||u- v||oo 

what finishes the proof. Here C, C and C" are positive constants. □ 

We have now all the keys for proving the existence of some symmetric invari- 
ant measure. Indeed we have just presented some continuous mapping which 
stabilizes a convex subset of the Banach space B. 

Theorem 4.5. There exists a symmetric invariant measure for l|l.ip . 

Proof. Let Mg defined by Lemma 14.31 Taking M > Mg, let us prove that 
A<^(Ca/) is a compact set. For this reason we shall estimate the following deriva- 
tive: 



, , 2 {F'*u)ix) + V\x) ^ 

u) (x) = — cxp 

e X^(u) 



{F'*u){y)dy + V{x) 



Let us analyze the different elements of this derivative. We have already seen 
in the proof of Lemma 14.31 that for any u e Cm the normalization factor Ac (u) 
satisfies 

-1- < C(6)^/M. (4.9) 

By (021), we obtain the bound: < J^' {F' * u) {y)dy < C2Mx^ (l + x^^). 
Furthermore by (V-4) and (V-7), we get some estimation oiV and its derivative: 

V{x) > Cix^ ~ C2X^ and \V' {x)\ < C, (1 + |xp9) for all a; e M. (4.10) 
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It remains to find some upper bound for the convolution term: \{F' * u) (.t) 
with X £ R+. By (F-4) and since u is an even function, 



\{F'*u) {x)\ = 



F' (x — z)u{z)dz 



< 



F'{x + z) + F'{x - z) u{z)dz 



J^\X + Z\ {1 + \X + Zp9-2) + \x-z\{l + \x~ z|29-2) ^u{z)dz 



Therefore: 



\{F'^u)ix)\ < [ C,22«-i||xp9-i + |z||x| 



2q~2 
2q-2\ 



By definition of ' 



-M, 



+ \X\ (1 + + \z\ (1 + |zp9-2) y[z)dz. 

for p > 4(j. Hence the moments of 



RI 



we have u{z) < x+j^ 
order 1, 2q — 2 and 2q — 1 are bounded: there exist some constants C and C", 
independent of the different parameters appearing in the system, such that 

* u) {x)\ <C{l + \x\ + + < C (l + . (4.11) 

To sum up: using (gJl), (|4.10p and (|4.1ip we obtain 



\{A'uy{x)\ < -C(e)VM(l + |x|29+i)exp 



{C^X^ - C2X^) 



Finally we deduce that there exists some constant Cc such that | (A^u)' {x) | < 
for all a; e M. 

Let us prove now that A'^Cm is compact. To this end, we take some sequence 
of functions (Mn)„gN '^m and focus our attention to the sequence {A'^Un)^^-^. 
According to the definition of A% for all x real the set {A'^u„(x),n S N} is 
compact. Furthermore the bound of |(A'^ti)' {x)\ is independent of the variables 
X and u € Cm- the equicontinuity condition for the appHcation of Ascoli's 
theorem is satisfied. Hence, we deduce that there exists some subsequence of 
A'^Un which converges to a limit function i; belonging to A'^Ca/. 
By Lemma [4.31 and Lemma [4.41 we can apply Schauder's theorem (Proposition 
14. ip for the operator A*^ on the function space Cm with M > Mq. We deduce the 
existence of some fixed point which is, by construction, a symmetric stationary 
measure for the diffusion Ijl.ip . □ 



4.2 Example: F{x) = f a;^ + f a;^ 

We have just shown the existence of a symmetric invariant measure for general 
self-stabilizing diffusions using fixed point arguments. Now let us study some 
particular case by a completely different way: the procedure shall be close to 
that developed in section [331 Let ^ be a potential satisfying (V-l)-(V-7). 
Let Ue be a symmetric invariant measure (Theorem 14. 5p . We denote by m2(e) 
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its second moment. The couple (m2(e), u^) is solution to some system like l|3.ip 
lla!2l). Indeed 



F * uJx) 



F{x — z)u^{z)dz 



2^ 



/? 4 3/3m2(e) 



with /3 > since F' is a convex function on IR+ . 

The expression dehmited by the brackets is just a constant so we obtain the 
following system of equations for m2(e) and u^: m2{e) = j^x:^v{m2{e),x)dx 
and Ue(x) = 1^(7712(6), x) where 



v{m, x) 



Therefore we introduce the function Xe (to) = x^h'{m,x)dx~m. By Theorem 
4.5\ we know that Xe admits at least one zero on R+ . Computing the derivative 
of Xe, we prove that the considered function is decreasing: 



exp 


-I (vix)+F{x) + ^x^^ 




rOO 

Jo exp 


>! 


{v{z) + F{z) + ^z^^ 


dz 



3/3 



x^vlm, x)dx 



x'^v{m, x)dx 



- 1< 0. 



The conclusion is immediate: there is a unique symmetric invariant measure. 
Obviously this result and the kind of method used to prove it are particular to 
our simple example. If the degree of the interaction function is strictly larger 
than 4 then it isn't enough to know the second moment in order to define the 
invariant measure: we need more moments and the proof of the uniqueness 
becomes awkward. 



4.3 Outlying invariant measures 

This section is essentially motivated by the uniqueness question for invariant 
measures. The existence of some symmetric measure was just proved in Sec- 
tion [4]TJ It suffices now to point out asymmetric stationary measures for self- 
stabilizing diffusions. In the general setting, the interaction function is polyno- 
mial: set F{x) =Y.l^,^x^\ 

Let u be the density of some probability measure with respect to the Lebesgue 
measure and /ii, • • • , fJ.2n-i denote its moments of orders 1 to 2n — 1 respectively. 
We assume they are finite. Then the difference D{x) := F * u{x) — F * u{Q) 
satisfies 

Dix) = F(x - a) - F{-a) + ^(a^p " E W^^^"''" 

■1—1 

= F{x -a)- F{a) + ^ ^^(Mp - a") [f^^Hx) - f'^p\o)) . 
p=i P- 
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Hence D{x) = Z„i{x) — Z„i{0) where 

2n-l 

I 

P 



2n-l , ,p 

Zmix) =Fix-a)+J2 -^(rnp ~ anF<^P\x). (4.12) 
p=i P- 

Since the convolution product can be expressed as a polynomial function which 
coefficients just depend on the moments of u, then the exponential expres- 
sion of invariant measure l|2.3p can be specified. Indeed equation (|2.3|) can be 
transformed into some system of equations whose unknown factors are the mo- 
ments of the measure. In order to introduce this system, let us define, for all 
fc € [1; 2ri — 1], the function 



ip^ [mi,-- - ,m2„-i) 



4 x'' exp [-f {V{x) + Zmjx) - ZmjO))] dx 

/^exp [-f (Vix) + Z^[x) - Zram\ dx 
J^x'' exp [-|W"„(x)] dx 



(4.13) 



/^exp [-|Wm(a;)] dx 

with the potential Wmix) = V{x) + Zm.(x). We construct the mapping: 

$(^) = (^(^\...,^i^\...,^(l,). (4.14) 

The measure associated to the density function u is invariant if and only if its 
moments vector (/ii, • • • ,/i2„_i) is a fixed point of the map 
We are going to show the existence of an asymmetric invariant measure defined 
by 2r7, — 1 parameters close to a, • • • , a^"~-^ respectively, in other words the 
outlying measure is close to the Dirac mass in the point a. More precisely, we 
shall prove that there exists a parallelepiped stable by which converges 

to the point (a, a?, - ■ ■ , a^"^^) as e tends to 0. As in the linear case, we shall 
proceed by applying the mean value theorem in order to obtain asymptotic 
developments in the small noise limit. 

Theorem 4.6. Let {i]e)e some sequence satisfying lim — and lime/rye = 0. 

£—►0 e—*0 

Under the condition 

2n-2 \p{p+2)( \\ 

J ^ a^<a + V"{a), (4.15) 

p=0 P' 

for any p > 0, there are at least two outlying measures uf and s 
for e small enough 

x''uf{x)dx - {±a)'' < p T],. (4.16) 
Proof. Let A > 0. Let us define the parallelepiped 

2ri-l 

C{e) = n y-paP'^Xrj,,aP+paP-'Xr,,]. 

p=i 
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Let m be an element of C(e) then there exist some coordinates (rp)i<p<2n-i 
which determine m through the equations nip = aP + rpT]^. By (|4.12|) and (|4.13p . 



we get 



/r 



-i(V{x)+F{x-a)) 



exp 



P— 



2n-l (-l)Pr-j 



dx 



dx 



We apply Lemma fA.TI and Remark lA.SI to the functions U{x) ~ V{x) + F{x — a)^ 
fix) =x'', fip = rp and Gp{x) = ^^F'^p^x). We obtain: 



ifl (m) 



ka 



fc-i 



'a + y"(a) ^ 



uniformly with respect to the coordinates {rp)p. By definition of the paral- 
lelepiped C(e) the coordinates satisfy \rp\ < paP~^X. Therefore, under condition 

633, 



ka 



fc-i 



a + V"{a) 



2n-l 



F(P+i'(a)| 



-pa 



p-i 



Since this estimate is uniform with respect to the coordinates, as soon as e is 
small enough, we have \ip'"f^\m) — a''\ < ka'^~^Xr]^, that means that ^^'^^(to) g 
C{e). 

Let us note that C (e) is a convex, closed and bounded subset of M^" ^ . Since 
the space dimension is finite, the continuity of ^^'^^ implies that the closure of 
the parallelepiped's image is a compact set. 

We can apply Schauder's Theorem (Proposition l4.ip and obtain that there exists 
some fixed point in the compact. In other words there exists m € C(e) such 
that the measure associated to the density 



exp [-IWM] 
/rGxP [-fW^m(^)] dz 



(4.17) 



is invariant. In a similar way, the measure defined by m~ is also invariant; here 
m^{k) = (— l)'^mfc. To conclude: we have at least two outlying measures, one 
around a and the second one around ~a. □ 

We can not prove at this stage the uniqueness of the couple of outlying in- 
variant measures (this question shall be explored in a subsequent work) . We can 
effectively imagine that other outlying measures could exist around a, around —a 
or even around other areas. Nevertheless we can develop a sharper description 
of one particular outlying measure: the measure close to Sa where S represents 
the Dirac measure. To do this it suffices to estimate its different moments, that 
requires the following preliminary result. 
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Lemma 4.7. There exists a unique solution (t^,-- - ,T2„_i) to the fo 
Cramer's system 

^ p! ^"^"^P^ fca^-i ' 4(a + y"(a)) 4a' ^ 
p— 1 \ \ // 

for 1 < k <2n — 1. This solution is given by 

T^" = fca , , / ,,,, \ l<k<2n^ 1. 4.19 

Proof. Let us denote by /2n-i the unit matrix of dimension 2n — 1 and for 
A e M^"~^, A""" represents the transpose of the vector A. Moreover we adopt 
the following notation (a;fc)i<fc<2n-i = {xi, . . . ,X2n-i)- The system l|4.18p can 
be written in this way: we define T = ('nc)i<fc<2n-i then 

(a + V"{a))hn-i + C,Cj] T = f ka'^-' ( ^ - ^) ) ^ 

. ^ ^' V \Aa + AV"{a)) 4a J J i<k<2n-i 



with the vectors C?" = {ka''-^)i<k<2n-i and Cj = (^^FC^+i) (a)) 
We define therefore 



l</£<2n-l 



^ (a + y"(a))/2„_i + CiCj. (4.20) 
Let us note that CiCjCiCj = {CjCi)CiCj and 



Since F"(0) = a, we obtain 



p=l ^ p=0 ^ 



A' = {a + V"{a)fhn^i + (2(a + V" {a)) + CjCi) C^Cj 

^{a + V"{a)fhn-i + (2(a + V"{a)) - F"(0))CiCj 

= (a + l/"(a))2/2„_i + (a + 2V^"(a))CiCj 

= (a + 2F"(a))A - V"{a) {a + F"(a)) /2„-i, 

We deduce that A is invertible, that is (|4.18p is a Cramer's system, and using 
l|4.20p we get explicitly the inverse: 

- yna)(aVy-(a)) ((- + ^^"("»^-^'-^^ 
^ (y"{a)l2n-i~CiCj). 



V"{a){a + V"{a)) 
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Therefore the Cramer's system i|4.18p admits a unique solution given by 



1 



V"{a){a + V"{a)) 

2n~\ 



^ p\ ^ ' Aa{a + V"{a)) 



ka 



fe-i 



4:aV"{a){a + V"{a)y 



V 



(a + y"(a)) 



(fc-l)F"(a)- J2 



p=2 



^^,_,ay(3)(a)-(fc-l)T/"(a) 



4a^^"(a)(a + \^"(a)) 



Indeed, we use 



(-1)^ 



■F(f+i)(a)aP-i = -F"(0) = -a, 



and 



2n-l 

E 

p=2 



(-IF 



F<'P+^\a)aP-^ =aF(3)(0) = 0. 



□ 



Theorem 14.61 points out the existence of two outlying measures, one con- 
centrated around a and an other around —a. According to Lemma [4.71 we get 
some sharper upper bound for the distance between Sa and some asymmetric 
invariant measure. 

Theorem 4.8. Under the condition l|4.15p . for any S > 0, there exists eg such 
that <f>'''^-' admits two fixed points with 

77i±(e) - ((±l)'=a'^ - (±l)V^"e) <S\T^\e, l<k<2n-l, e< eo- (4.21) 



Proof. It is similar to the proof of Theorem 14.61 

Let (5 > and C(e) = Ul^i^K ~ i^p +paP~^d)e,aP - (t^ - paP-^S)e]. We 
choose an element m in the parallelepiped C(e). For all 1 < p < 2ri — 1, there 
exists a coordinate Sp G [—^',5] such that nip — qP — (t° + paP~^5p)e. By l|4.12p 
and (|4.13p . we obtain 



Ir exp 


'2ET=i'^ir'p+paP-%)F(P)i.) 






2Ep:T' ^-^{r^p+paP~^5p)FiP){x) 
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We apply Lemma IA.5I and Remark IA.8I with the following functions: U{x) = 
V{x) + F(x - a), = rO +paP-^5p, G = and fp{x) = 2^^F^p'>{x). Hence 



4(a + y"(a))2 
(-IF, 



ay(3)(a) _ (a + V^"(a))((fc- 1) 



-4a J2 ^-^{r'p+paP-'Sp)F(P+'Ha))\e + oie) 
p=i P' 



1 



-n 



a + y"(o) L4(a + y"(a)) 



k{k- l)a 



fc-2 



This estimate is uniform with respect to the variables {Sp)p. 

We denote by d| the difference \(pl^\m)-~a'' + T^e\. We compute this expression: 



4< 



4(a + V""(a))2 a + y"(a) ^ p: 

p— 1 

fc(fc - l)a'--2 
4(a + y"(a)) 



-1 ^""-^ (- 



—5:^ V ^~ ' S F'^P+'^'>(a)aP-^ 

a + V"{a) ^ {p-iy " ^ ' 



e + o{e). 



According to the Lemma [4.71 and using the condition (|4.15p . we obtain, for e 
small enough, 



< S 



a + V''(a) ^ p\ 



In other words, $(^^(m) e C(e) in the small noise limit. The application of 
Schauder's Theorem (Proposition I4.ip permits to prove the existence of some 
fixed point in the compact. Therefore there exists m € C{e) such that the 
associated measure Ue^m(x) defined by l|4.17p is invariant. In the same way, 
the measure defined by m~ is invariant with m~{k) = (— l)'^mfe. Finally the 
continuous map admits two fixed points m'^{e) satisfying l|4.2ip . □ 

Remark 4.9. 1. In the particular case: F'^P^a) > for allp € N, the condition 
for the existence of outlying measures becomes V" {a) > Fg'(2a) where Fq is 
defined by F{x) = f + Fo{x). 

2. In the linear interaction case: F{x) = ^x"^ , (|4.15p is equivalent to the simple 
condition V"{a) > which is in fact always satisfied according to (V-3). In 
other words we obtain the existence result presented in the linear interaction 
case. 
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A Annexe 



We shall present here some useful asymptotic results which are close to the 
classical Laplace's method. 



Lemma A.l. Let M > 0. Let us assume that U is C'^{[M,oo[)- continuous, 

(U'(x)) 



U{x) 7^ and U"{x) > for all x G [A/, cx)[ and limj;_>oo rn'M)^ ^ ^- V 



e ^^^'^ is integrable on R then for any m e R; 

+00 p-U(x) px U(x) 

e-^^'Ut « and / e^^'Ut « |— - as x ^ oo. (A.l) 



X 



U'ix) J„, U'{x) 



Proof. Since x e"'^^^' is integrable and since these properties are satisfied: 
U{x) ^ and U"{x) > for x > M, we know that lhnx^ocU{x) = +oo. 
Furthermore there exists some Mq > M such that U'{x) > for a; > Mq. Hence 
for t > Mo we obtain 



Therefore 



u'{t)) (c/'(t)r' 



U'ix) J, {U'{t)f 
Using the assumptions of the statement we have (u'U))'^ e'^'^^^dt > 0. Hence 



X {u'{t)y 

iu'{x)r 



I < e ^^^^U' [x) ^ . Moreover lim^r ^±00 (^'(i^ja =0- As a consequence for any 

(5 > 0, there exists Mi{5) > Mq such that (1 + S)I > ^ipj^- The estimation 
of / can be deduced easily. The second equivalence can be obtained by similar 
arguments. □ 

Lemma A. 2. Set e > 0. Let U and G two (M) -continuous functions. We 
define ~ U + f-iG for f-i belonging to some compact interval X of R. Let us 
introduce some interval [a, b] satisfying: U'^{a) ^ 0, U'^{b) 7^ and U^{x) admits 
some unique global minimum on the interval [a, b] reached at x^ s]a, b[ for all 
^ <E T. We assume that there exists some exponent kg independent of ^ <E 2 
such that 2ko = min,eN* ^uj^^x^j^) ^ o|. Then taking the limit e we get 



-(l + ox(l)), (A.2) 




where T represents the Euler function and oi(l) converges towards uniformly 
with respect to ^ 
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Proof. We define 77^ = ''^^fcp)! ^ • Let us note that 77^ depends continuously on 
fj,. Since is regular and admits some unique global minimum for x = x^ , there 



exists tq > independent of the parameter ^ such that tq < min {x^ 



for all ^ e X and such that the minimum on the interval [a; — r] [J[x^ + t; b] 



denoted by iZ^(T) is reached on the boundary { 
Consequently 



exp 



dt 



cxp 



Defining U = J^"_r ^xp 



Xa — r; + r} for all r < tq. 



dt < (b — a) cxp 



1^ - It\ < {b- a) cxp 



dt, we obtain the following bound: 



(A.3) 



Let us first estimate Ir- By the mean value theorem, there exists some constant 
C > independent oi ^, <E 2 such that, in a neighborhood of x^, the following 
bound is satisfied: |J7^(t) - C/^(a;^) - 7?^(t - x^J^fco | < C|i - Hence 



Ji exp 



< — exp — — — < Ji exp 

^ 6 



(A.4) 



where 



Jr = 



cxp 



2fco Jo 



by the change of variable t := (^^) {t')'^>'o . A simple integration leads to 



< 



.21-0 



tw-ig-'di - r ( ) < 0. (A. 5) 



In order to conclude we choose a particular value for the variable r namely 



exp 



log(e) 



Then we get: for C G K, Z > 0, 



lim e 

€^0 



^2fco + l 



1, lime ''f = — i = lim e e 

£-►0 e2fc5-"^ 



0. 



These convergences are uniform with respect to the parameter Applying these 
asymptotic results to ijA.Sp . (|A.4| and ijA.SP permits to prove the statement of 
the lemma. □ 

Lemma A. 3. Let U and G he two C°°([a, b]) -functions. We define = U 
for fj. belonging to some compact interval J of M. We assume that admits a 
unique global minimum on the interval ]a; b[ reached atx = Xfj_, with U'^{Xfj_) > 0. 
Let fm be a -continuous function for any parameter value m belonging to 
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some compact set A4. We assume that there exists some constant A such that 
\fm\x)\ < A for all m G A4, x G [a,b] and < i < 3. Then the following 
asymptotic result holds: 



fm(t)e - dt = . — e 

V 1^2 



|/m(%) +7o(M)e + oiA^(e)| (A.6) 



with 



loilj) = fm{x^) 



48 



16 Ul 



Us 



4U2 



(A.7) 



HereUk - ui''\x^) and ox»(e)/e converges to Q as e becomes small uniformly 
with respect to the parameters m and ^. 



Proof. First we split the integral into two parts: 



fm{t)e " dt 



fm{t)e ' dt 



h+h 



with some arbitrary p > which should be specified in the following. 

2 

Step 1. We shall prove that the second integral is negligible as ^ - 



means that I2 — oxm{^^^^s 



}. We get 



l2<{b-a) sup |/m(z)|exp 



-2- 



inf- 



llx^-p-.x^+p]" U^{z) 



00 that 



(A., 



Since the global minimum of C/^ is unique and due to the regularity of Uf^ with 
respect to the parameter fi, we deduce that the minimum of the function on the 
interval [a;^ — + p]'^p|[a;6] is reached on the boundary provided that p is 
small enough. The development Ufj,{x^ ± p) = [/^(x^) + \U'i[{Xfj)p'^ + ox{p^) 

implies, as already claimed that I2 = ojm as p^/e 00. 

Step 2 . Let us focus our attention to the integral on the domain [a;^ — p; + p] . 
The function fm can be developed in the neighborhood of : 

fniix) ^ f,niXf,) + f!^{Xf,){x~Xf,)+^f!^{Xf,){x~Xf, f + ^fj'^\w,nA^))i^-^f^f 

with the value Wm,p{x) between x^ and x. Taking into account these different 
terms, the integral Ii can be split into 4 different integrals respectively lo,...,!^. 
For each integral we shall analyze the asymptotic behavior. 
Step 2.1. Asymptotic behavior of I3. By definition Wm^^{t) e [a;^ — p; + p] 

when t £ [Xfi — p;Xfi + p]. Moreover, by assumption Ifm {wm^f_i{t))\ is upper 
bounded by some constant A > independent of m and p. By Lemma IA.2I 
applied to 2f7^ , for p < 1 and e small, we obtain the existence of some constant 
C > 0, independent of the parameters m and p, such that 



1 < ^p' 



{x^+l)Ab 



(3;^ — l)Va 



e ~^~dt < Cy/np^ 



U''{x^) 
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Hence, if p = o(e) then the following asymptotic result holds 



(A.9) 



Step 2.2. Asymptotic behavior of /2. Using the C'^-regularity of C/^ that is 

U^it) = U^ix^) + ^U';{x^){t-x^)^ + ^ui^\y^it)){t-x,,f with y^it) belonging 
to [x^ - p;x^+ p], we get 



-'2 — e 



{t - Xf_,) e 



(3) 

Since belongs to some compact set, the third derivative C/^^ (^^(i)) is 

bounded by some constant independent of p. Applying the following change of 
variable u ~ {t — Xp)^C/^'(x^)/e yields 



J2e 



< 



2/2 e 



with J2 = /o^'"^''^-''^ y^e-^'du. If 4 and ^ 



< J2e 



00 then 



I■2 = ^/^ 



(A.IO) 



Step 2.3. Asymptotic behavior of Ii. Let us develop the function C/^ in 
the neighborhood of x^: C/^(t + x^) = Uf,{xo) + ^U'/^{x^)e + \u'i?\x^)t^ + 

^U'lf' {y^{t))t'^ where S [x^ - p, x^, + p]. The regularity of U^{x) with re- 
spect to both X and /i implies the existence of some constant C > independent 
of p which bounds the forth derivative of on the integral support. Therefore 
we have 



with Jp = ze 



dz and Uk = Uli (xp). Since |e ^ ~ 1 + x 



< 



xpe'^l, we deduce that, for any z € [— 



e-^-^ - 1 + ^ 
3e 



3 Ulz^ 



18e2 







< 






3 



3 Q o 
P« |»3lp3 

; p 3e 

£3 



We define mp{l) — z'e ' ^ dz and np(^) = j^^ |z|'e dz. Some estima- 



"2 ^2 



tion of the integral Jp points out directly: 



Jp - mp(l) + -^"ip(4) 



18e2 



Jo 




< 2 








3 



3 Q o 

-36 3. np(l). 
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Symmetry arguments permit easily to deduce that rnp{l) = mp(7) = 0. Finally 
it suffices to compute rnp{4) and np(l). To this end we introduce the change of 
variable u := and let p^/e tend to infinity: 



3x/^ / 1 



4 \u';ix,) 

IS 9 3 

To sum up: if ^7 ^ (that is |j- = o{e2}) then 



£2(1 + 01(1)) and np(l) = -——(1 + 01(1)) 



U;iixp)J ^^e-^(l + o.(l)). (A.ll) 
Step 2.4. Asymptotic behavior of /q- Let us first study the following integral 



cxp 



3e 



12e 



We recall the usual notations Uk ~ u'lj^\xp). The arguments are similar to 

2 

those used in Step 2.3. Since e^" - 1 + u - ^ 
we get 



< |upel"l, for any z e [-p; p] 



e T2r^ - 1 + —z^ + — z - 



3e 12e 2 V 3e 12e 



-z" + ——z ' 



Adopting the same notations as in Step 2.3 and using symmetry properties, the 
following bound (uniform with respect to the parameter p) yields 



4--p(0) + ^-p(4)-^(^^) -p(6)-^(^) m,(8) 



1 fU: 



12e 



2 V 3e 



1 I Ua 



2 V 12e 



By the usual change of variable u := iji^iiiil^;^ -^ve emphasize some asymptotic 
estimation of /q as p^/e co and p'^/e ^ 0: 



V C/;(.T^) \ 16C/;(x^)2'^ 48C/;(x,)3' + °^^'^ 



" ^'C/;(.T^)[^ 16C/;(x^)2- 48C/;(x^)3' 
We apply the mean value theorem to the function Uf^: 

Up{x, + z) = Up{x,) + + ^z^ + + ±-^U^^)iy^it))z\ 

with y^(t) S [Xf_i — p,Xfi + p] and |z| < p. From this equality we deduce an 
estimation of the distance between the integrals Iq and Iq. 
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We denote hy V = e ' Iq — fmixfj)!^ this distance. Then there exists some 
constant C > independent of yU and m such that 

rP 



PI < \Im{x^)\ f e 



"2 2 M3 3 »4 4 



1 _e-6fcf^'"Cyf.(^+^M)^' 



- 60e P 1/ 
If both conditions p^/e ^ 00 and p^/e ^ are satisfied then the integral term 



o e 



in the preceding inequahty is obviously equivalent to ^xpr^r^- The following 
equivalence holds for the initial integral /q: under the assumption that ^ = 
we get = oxM (^^) ^^'^ consequently 



V 1 16 Z^l 48 Ui J 

Step 3. To sum up: in Step 1, we proved that it suffices to estimate the 
integral Ii which can be split into 4 terms. Each of them has been estimated 
in equations (|A.9p . (|A.10p . (|A.lip and l|A.12p . The whole integral has the 
asymptotic equivalence (jA.6P as soon as p^/e — > 0, p^^/e'^ — > and p^/e^ — *■ 0. 

g 

The particular choice p = e^o fulfills all these conditions. □ 

We can extend the statement of the preceding lemma to integrals with un- 
bounded supports. 

Lemma A. 4. Let U and G he two C°° (M.)- continuous functions. We define 
Ufj, = U + p.G for the parameter /i belonging to some compact interval X 0/ R 
We assume that U^{t) > for \t\ larger than some R independent of p and 
that Ufj, admits a unique global minimum at with [/^'(x^) > 0. Let fm be 
a -continuous function depending on some parameter m which belongs to a 
compact set A4. Furthermore we assume that there exists some constant A > 
such that |./m(<)l < exp [A|J7p(t)|] for all t > R, p e J, m e M and \fm\ is 
locally bounded uniformly with respect to the parameter m e Ml for < i < 3. 
Then the following asymptotic result holds as e tends to 0; 

/m(Oe dt = e Y ^ + 7o(Ai)e + oiAi(e)|, (A.13) 

where 70 (/i) is defined by (|A.7P and oim (e) / e converges io as e — > uniformly 
with respect to the parameters m and p. 

Proof. Let R > such that U^{t) > for t > R. The initial integral can 
be split into two integrals: the first one denoted by Ii concerns the compact 
support [— i?, R\ and the other one I2 concerns the complementary support. For 
/i it suffices to apply Lemma [A. 31 in order to get the asymptotic development. 
It remains then to prove that I2 is negligible with respect to Ii that is I2 = 
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£2 6 ; >. Using the change of variable t := (- — A) ^ s the following 



bound holds: 

r+oo 

Ihl < 2 / exp 



2 - Ae 



+C30 

21 



c^i < 2W- — / exp -s^ 



Lemma Fa. II permits to prove as claimed that I2 can be neglected. □ 

Lemma IA.4I can be applied to particular functions fm ■ 

Lemma A. 5. Let U and G he two C°^{R.)- continuous functions. We define 
Ufj_ = U + i-iG with /i belonging to some compact interval T of M. We as- 
sume that U^{t) > for \t\ larger than some R independent of fi and that 
Ufj, admits a unique global minimum at with U'^{xp,) > 0. Let fm be a C^- 
continuous function depending on some parameter m which belongs to a compact 
set A4 . Furthermore we assume that there exists some constant A > such that 
\fm{t)\ < X\U^{t)\ for all t > R, j.1 e I, m G M and that \fm\ is locally 
bounded uniformly with respect to m <E Jv[ for < « < 3. Then, for any n > 1 
and asymptotically as e we obtain the estimate 



J^t-ef-('h^^dt 



nx,, 

X?. ^ 



n-2 r 



J^ef^i^)e^dt 



^ - ?i + 1 - 2xf,f^{x^,) 



oiAi(e), 



where Ui = uIjI\x^) and oxm{^)I^ converges to as e ^ uniformly with 



respect to the parameters m and ^ 



Proof. We just apply two times Lemma [A. 41 the first time to the denominator 
D'^ that is for the function t e^-mit) ^^^^ ^j^g second time to the numerator N'^ 
for the function t t'^e^"^^*\ The following asymptotic result holds 

= e ^ J— e^'"(^")|l +7de + oiM(e)| (A.14) 



where 



The numerator normalized by x^ i.e. /x^ satisfies some similar identity as 
D'^ , namely (|A.14p with 7^ replaced by 7„: 

^" " l6Ui J {x^ ^ /rnl^M^J 4^2 

+ (^^^^ + 2^f^{x,) + + f:jx,)A 

\ xf, X^, J AU2 

The estimation of the ratio is then a classical exercise of asymptotic analysis. □ 
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The next lemmas are generalizations of Lemma IA.4I and Lemma IA.5I to 
functions G depending on the small parameter e. 

LemmaA.6. LetU and G be two C°°{R)- continuous functions such that U{t) > 
t^ for \t\ large enough and \G{t)\ < X\U(t)\ + C for some constants A > 
and C > 0. Moreover we assume that U admits some unique global minimum 
reached at xq with U"{xq) > 0. For any sequence (r]e)e satisfying lim^-to f?e = 
and linie^o ^/Ve = we define U^^^ = U + rj^jiG depending on the parameter ^ 
which belongs to some compact interval 2 o/M. Let f a -continuous function 
such that \ f(t)\ < e^l^*^*^! for all \t\ large enough and such that \ fm\ is locally 
bounded uniformly with respect to m ^ Ai for < z < 3. Then, there exists 
eo > such that the potential Ue,^ admits a unique global minimum reached at 
Xe.i_i for all e < eg. Furthermore the following asymptotic results hold 



^U"ixo)- 



=xo- M77777rT^<^ + 02:(f?£) (A. 15) 



f{t)e —dt^J——-e : {fixo)+j^Ve+oxive)), (A.16) 

I y U (Xo) ^ ' 

where 

and ox{ri()lr\^ tends to as e ^ uniformly with respect to the parameter fi. 

Proof. Let us first prove that the potential C/c^^(a;) admits a unique minimum 
for X = Xe^fj, with limj-^o^e^^ = Xq. By the definitions of {r]e)e and U^^fj,, the 
following convergence holds 

limU,^f,{xo) = U{xo). (A.17) 

Since xq is the unique global minimum of U, for any small R > there exists 
pR> such that inf^g[a;(,_/{_2,„+^]c U{x) > U{xo)+pR. We deduce the existence 
of two small constants and eg such that 

U,,^{x) > (1 - fiX7j,)U{x) ~ r],fiC > U{xo) + p'r, (A.18) 

for all e < eg and x £ [xq — R,xo + R]'^. By ijA.lTP and ijA.lSP we obtain: for 
any R> the global minimum of the parametrized potential C/e,^ is reached in 
the interval x £ [xq — R, xq + R] provided that e is small enough (uniformly with 
respect to p). Moreover this global mimimum is unique. Indeed U"{xq) > and 
the regularity of U implies that U" (x) > for all x in some small neighborhood 
of Xq. Since C/"^ converges towards U" as e uniformly on each compact 
subset of R, we obtain that J7"^ (a;) > for all x £ [xq ~ R, xq + R] provided 
that R and e are small enough. The minimum is actually unique, we denote its 
locaHzation x^^^ and point out that, for e small, U'^_^{Xf_^^) > uniformly with 
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respect to yU. 

Let us determine x^^^. By applying the mean value theorem to Ue.fj., we get 

= U'^,^{Xe,f,) = U'{xo) + fJ.%G'{xo) + U"^{x){Xe,f_t - Xo), 

where x is in between a;o and x^^fj,- Since the second derivative is continuous, 
i7"^(i) is uniformly bounded. Moreover U'{xo) = 0. Consequently Xe,^ ~ xq = 
Ox{r]e). Using the same argument for the second order asymptotic development 
of U'^^^{xe,fj.), that is 



o = u- 



'(xo)+WeG'(xo) + (t/''(xo)+A^r/eG''(xo))(z,,^-xo) + ^^^|^(xe,,.-xo)^ 

we obtain the announced estimate (jA.lSp . Finally let us prove the estimate 
IjA.lCP . The statement of Lemma [A.4l can be applied to C/e^^ since the asymptotic 
result IjA.lsp is uniform with respect to the parameter /z. So it suffices to 
consider the case when fi is replaced by 11%. We immediately obtain 



f{t)e-^dt ^ ^) /(xe,^)e"^^^^ (1 + oxfa)) . (A.19) 



f{x,^^) = /(.To) - mJ\xo)jyj^Y^ + oi{ri,), 



It remains to approximate f{xe,^i) and C^"^(xe^^) using (|A.15|) . Due to the 
regularity of both / and U , the following developments hold 

G'(xo) 
'U"{x^) 

U'l^ix,^,) = U"{x,) + ^,r^,(G"{xo) - U^'\^o)§^) +oi(r?.). 

The statement of Lemma [A.6I is obtained just by combination of the two pre- 
ceding asymptotics and (|A.19|) . □ 

We are now able to present a statement similar to Lemma IA.5I for some 
potential C/^ depending on the small parameter e. It suffices to consider a ratio 
of two integral terms. Then an immediate application of Lemma [A.6I leads to 
the following result. 

LemmaA.7. LetU and G be two C°°{R)- continuous functions such that U{t) > 
t^ for \t\ large enough and \G{t)\ < X\U{t)\ + G for some constants A > 
and G > 0. Moreover we assume that U admits some unique global minimum 
reached at xq with U"{xo) > 0. For any sequence {q^)^ satisfying lim^ ?7c = 
and lime_»o e/?7e = we define Ue,^, = U + rj^iiG depending on the parameter ^ 
which belongs to some compact interval! ofM.. Let f a -continuous function 
such that \ f{t)\ < e^l^(*)l for all \t\ large enough and such that \ fm\ is locally 
bounded uniformly with respect to m e A4 for < i < 3. Then as e ^ 0, we 
obtain the following estimate 

fixo)G'ixo) 

where ox{r]e)/r]'^ tends to as e ^ uniformly with respect to the parameter ^. 
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Remark A. 8. The statements of Lemmas \ A. ^A.7\ can be easely generalized, 
replacing the parametrized function = U + fiG by Ufj, = U + l^iGi where 
/Li = (/Lii, . . . , fXk) € Ii X . . . X Xfc. The convergence results are then uniform with 
respect to all parameters. 
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